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A new functional form for the exchange enhancement in the generahzed gradient 
approximation within density functional theory is given. The functional form satis- 
fies the constraints used to construct the Perdew-Burke-Ernzerhof (PBE) functional 
but can be systematically varied using one parameter. This gives the possibility 
to estimate the reliability of a computational result or to fit the parameter for a 
certain problem. Compared to other semi-empirical functionals, the present has the 
advantage that only one physically transparent parameter is used and that the fit- 
ted functional will obey the same exact conditions as PBE functional. Furthermore 
the simple form of the exchange enhancement means that oscillating terms in the 
exchange potential are avoided. 

PACS numbers: 71.15.Mb,61.50.Ah 

I. INTRODUCTION 

Density functional theory (DFT) has made it possible to calculate ground state properties 
of even very large systems efficiently and accurately. As the exact functional is unknown, 
practical applications of DFT require an approximate exchange correlation energy, Exc[n\.— 
Of these, the local spin density approximation (LDA)ii^i^ has been successful despite its 
simplicity. The generalized gradient approximations (GGAs) are an attempt to improve on 
the LDA. It is possible to construct a GGA free of empirical parameters, known as the Perdew 
Burke Ernzerhof (PBE)-GGA, by demanding that the GGA obeys certain fundamental 
constraints.- Numerical tests have shown that the PBE-GGA gives total energy dependent 
properties in good agreement with experiment .-^^ Consequently, the PBE-GGA functional 
has been extremely influential, both for performing actual calculations and as a basis for 
functionals involving higher derivatives and exact exchange.- 
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As both the density and the gradient can only be constant in the homogeneous electron gas 
limit, there can be no unique GGA- and the constraints of the PBE-GGA are not sufficient 
to uniquely define the functional. The PBE functional form was based on a numerical 
GGA^'^° where a model of the exchange correlation hole was constructed to satisfy known 
exact hole constraints. The constraints were satisfied using a sharp real space cutoff and 
a damping function, which were choices of the authors and different choices would lead to 
different functionals.— This was recently demonstrated by Wu and Cohen (WC),^^ who used 
the gradient expansion for slowly varying densities^^ to construct a GGA with a functional 
form corresponding to a diffuse cutoff.— The WC-GGA damps the gradient enhancement, 
which results in improved equilibrium volumes for densely packed solids, but poorer exchange 
energies of atoms compared to the PBE.'^i Another modification of the PBE, the RPBEj^, 
chooses a functional form with a larger gradient enhancement, which generally leads to very 
good atomic exchange energies but poor equilibrium volumes of solids.^ The PBE, RPBE 
and WC functionals differ only in the functional form of the exchange energy enhancement 
and all satisfy the same conditions as the PBE. 

A different approach to has been to construct GGAs as parameterized fits (See Kurth et 
al.- for an overview). These semiempirical functionals are often precise for a certain set of 
compounds, but the parameters are rarely physically transparent and it is consequently not 
clear when they could fail. Furthermore they generally violate the known exact conditions 
satisfied by the PBE-GGA. This paper introduces a new functional form that brings the 
WC, PBE and RPBE exchange functionals onto a common ground and allows semiempirical 
functionals, that do not violate the constraints of the PBE-GGA, to be constructed. The 
shape of the functional is controlled by one physically transparent parameter and in the 
construction of the functional it is ensured that the exact conditions of PBE are not violated. 



II. BACKGROUND 



The GGAs write the exchange energy density pr. particle as 

e^(n, s) = e^.(n)F^(s) (1) 

where e^in) is the LDA exchange energy density pr. particle. is the enhancement factor 
due to density gradients and is dependent on the reduced density gradient, s = — — ^YT^yrijs- 

2(37r ) ' Tifrj 
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The PBE enhancement factor is given as 

i^r^(^) = i + Mi-^^) (2) 



1 + x/ K ^ 

where x = fip and p = s^. The parameters of Eq. ([2]), /i = 0.2195 and k = 0.804, are 
determined to ensure that the exchange gradient correction cancels that for PBE correlation 
as s — i> and to ensure that the local Lieb-Oxford bound is obeyed.™ The RPBE functional 
differs in the functional form 



j^RPBE , 



X) 



1 + e""/") (3) 



while the WC functional retains the functional form of PBE but introduces a complicated 
^wc ^ io/81p+ (/i- 10/81)pe-P + ln(l + cp2).ii The enhancement factors of WC, PBE and 
RPBE are shown in Fig. [H 

Obviously the WC, PBE and PRBE functionals satisfy the local Lieb-Oxford bound. By 
Taylor expanding for p it is clear that all three exchange functionals have the correct 
first order expansion, F^ = 1 + fJ^p + O (p^), and thus obey the cancellation condition: 
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l+^p-^p^ + ^p^ + 0{p') (4) 

Fr'' = i+,p-l^p'+-^y+o{p^) (5) 

F^'^ = l + f.p-(^ + f.-b-c)p' 

^ ^ (6) 

+ K + ''^~''''''' ^'-i^)p'^oip^) 

\ /t^ Kj Zs I 

The exchange enhancement, which will be named PBEa, that will be presented here has 
the form 

:^PBEa ■ I ■ 1 



pPBE. ^ 1 + ^ M - , (7) 

where a is a parameter. The Taylor expansion 

i^f^^" = 1 + - + + + + o (/) (8) 

lot K bK'^a^ 

shows that PBEa satisfies the cancellation condition for all a. The PBE and RPBE are 
then special cases of PBEa = 1 and PBEa oo respectively. The PBEa functional is 
shown for two different a in Fig. [T] where it is seen how the size of the gradient enhancement 
can be controlled through the a parameter. Other attempts have been made at systemati- 
cally varying the PBE functional form either by making an expansion of the enhancement 
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FIG. 1: (color online) Exchange enhancements, F^, due to density gradients. The PBE, RPBE 
and WC functionals are shown with thick lines. Two PBEa functionals are shown with thin lines 
and marked with the a used. 

coefficients^ or by varying the k parameter.— However, these functionals either violate 
the local Lieb-Oxford bound or can only ensure it through a fitting constraint. 

III. DISCUSSION 

In the following, the influence of the parameter a will be discussed. All calculation have 
been performed using the WIEN2k code.— The basis sets and fc-meshes have been converged 
and equations of states (EOSs) have been found by fitting the stabilized jellium EOS''^ to 
the calculated total energies. 

A. Influence of a 

The influence of the parameter a can be illustrated by comparing the body centered cubic 
(BCC) phase of Fe with the hexagonal close packed (HCP). The example is well studied 
and it is known that LDA wrongly predicts the non-magnetic FCC phase to be stable at 
zero pressure.— EOS parameters have been calculated for BCC and HCP iron using LDA 
and PBEa for three choices of a, Table [H Using the EOS parameters, BCC-HCP phase 
transition pressures of 14.0, 16.4 and 18.3 GPa are predicted for PBEa = 0.8, a = 1.0 and 
a = 1.2 respectively, in reasonable agreement with a recent experimental study which fixed 
the transition pressure between 14 and 16 GPa.— 



Lieb-Oxford bound 
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P T^BCC 




^0 


^0 








(niRy/atom) 


(aQ/atom) 


(GPa) 




BCC 


LDA 




0.00 


70.19 


250.69 


4.97 




PBEa = 


0.8 


0.00 


76.18 


190.69 


5.29 




PBEa = 


1.0 


0.00 


76.45 


188.11 


5.32 




PBEa = 


1.2 


0.00 


77.03 


182.34 


5.33 


HCP 


LDA 




-9.76 


64.80 


337.72 


4.30 




PBEa = 


0.8 


6.15 


68.81 


288.06 


4.52 




PBEa = 


1.0 


7.25 


69.10 


285.63 


4.53 




PBEa = 


1.2 


8.10 


69.30 


282.70 


4.54 



TABLE L Equation of state parameters for BCC and HCP iron. 

The iron example shows how the a parameter can be interpreted. A small a favors the 
close packed structure (in the extreme case of the LDA the HCP structure is predicted to be 
stable at ambient pressure) while a large a favors an open structure. This also illustrates the 
dangers of the most obvious use of the PBEa functional: Fitting a for a certain problem. 
If one uses data set of densely packed structures and a small a, one risks biasing the results 
towards the densely packed structures. If one uses a data set of open structures (or in 
the extreme: atoms) and obtain a large a one risks biasing the results towards the open 
structures. Often it is probably better to use the well tested PBEa = 1.0 functional. 
However, varying a will provide information on how sensitive the results obtained are too 
variations of the functional. Such a philosophy was also recently proposed in a statistical 
study of the exchange enhancement.— The PBEa functional would probably better suited 
for this approach than the form actually used,— which was similar to the mPBE.— 

B. Limits of a 

Another question would be whether there is an upper and an lower bound for a. The 
good performance of the RPBE for atoms and WC for densely packed solids can give the 
impression that these functionals form the upper and lower bound of the gradient enhance- 
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X-OEP 


WC 


a = 0.52 


a = 1 


a = 2 


a = 5 


a = 20 


RPBE 


He 


-1.0258 


-0.9805 


-0.9916 


-1.0051 


-1.0145 


-1.0212 


-1.0249 


-1.0262 


Ne 


-12.1050 


-11.8676 


-11.9597 


-12.0275 


-12.0716 


-12.1015 


-12.1176 


-12.1231 


Ar 


-30.1747 


-29.6846 


-29.8646 


-29.9814 


-30.0551 


-30.1039 


-30.1299 


-30.1388 


Kr 


-93.8330 


-92.8559 


-93.1862 


-93.3769 


-93.4925 


-93.5671 


-93.6060 


-93.6193 


Xe 


-179.0635 


-177.5052 


-177.9762 


-178.2368 


-178.3923 


-178.4914 


-178.5428 


-178.5602 


MARE 




1.98 


1.37 


0.85 


0.50 


0.26 


0.18 


0.16 



TABLE II: Atomic exchange energies in Ha. The X-OEP marks the optimized effective potential 
values taken from Kurth et al.— MARE is the mean absolute relative error in %. 

ment. The RPBE performs extremely well for the exchange energy of atoms which would 
indicate that there is no upper limit for a. This is shown in Table HTl where MARE decreases 
with a. 

If exponentially decaying densities in atoms constitute one extreme, the other extreme is 
the slowly varying density. For the slowly varying density the expansion of is known^^ 

Sr 2025^ 405^-^ \ > \ ) 

where D ^ and q = = ^^g^^^^z/" „5/3 • The expansion in Eq. ([9]) was the inspi- 

ration for the WC functional,— where the Laplacian terms were avoided by introducing 
the approximation: q 2/3p. Introducing this approximation into Eq. ([9]), one obtains 
pSvB _ ^ ^ l0/81p - 0.0881207p2. From Eq. (E]) it is seen that if one sets a = 0.52 in 
the PBEa functional the same second order term is obtained. In Tables IIIII and IIVI the 
equilibrium volumes and bulk moduli of 18 solids^*^"^ calculated with the PBEa = 0.52 
functional are given and compared to the results of the LDA, PEE, RPBE and WC func- 
tionals. It can be seen that the PBEa = 0.52 functional gives results that are similar to the 
WC functional, being slightly worse on the equilibrium volumes and slightly better on the 
bulk moduli. One cannot exactly reproduce the WC functional with the PBEa because of 
its different X2 but both the WC and the PBEa = 0.52 functional were constructed to give a 
good description of the known gradient expansion of for a slowly varying density.— The 
performance of PBEa = 0.52 thus confirms that this expansion is important for densely 
packed solids.— 
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Solid 


u 






yRPBE 




yPBEa=0.52 


Li 


141.83 


128.32 


136.68 


142.10 


138.37 


134.13 


Na 


255.48 


224.17 


249.50 


266.49 


249.48 


241.37 


K 


481.31 


432.57 


501.21 


545.80 


497.04 


479.00 


Al 


110.59 


106.66 


111.37 


113.50 


109.66 


110.16 


C 


76.57 


74.50 


77.03 


78.13 


75.89 


76.36 


Si 


270.11 


266.22 


274.41 


281.58 


270.79 


273.59 


SiC 


139.64 


136.93 


141.88 


144.06 


139.56 


140.58 


Ge 


304.61 


300.20 


322.97 


332.33 


309.20 


316.80 


GaAs 


303.96 


297.43 


320.53 


330.42 


306.63 


314.20 


NaCl 


295.48 


275.22 


312.64 


337.12 


300.86 


300.10 


NaF 


165.18 


153.52 


175.11 


188.37 


169.36 


168.36 


LiCl 


224.58 


207.45 


231.84 


247.88 


220.74 


223.53 


LiF 


108.79 


100.70 


113.23 


120.31 


108.72 


109.41 


MgO 


125.62 


121.69 


130.35 


134.72 


126.71 


127.93 


Cu 


78.91 


73.56 


80.58 


83.44 


76.74 


78.80 


Rh 


92.43 


89.95 


95.29 


96.94 


92.45 


94.11 


Pd 


98.62 


95.82 


103.19 


105.93 


99.05 


101.43 


Ag 


113.66 


107.88 


119.31 


124.57 


112.49 


116.19 


MRE 




-5.26 


2.88 


7.32 


-0.03 


0.52 


STD 




3.30 


2.76 


4.22 


1.70 


2.56 


MARE 




5.26 


3.54 


7.32 


1.32 


1.93 



TABLE III: Volume of the primitive unit cell (in ag). MRE is mean relative error. STD is the 
standard deviation of the relative error. MARE as in Table HIl All errors are in % 

C. Potential 

Though often ignored in the construction of a functional, the resulting exchange- 
correlation potential is an important quantity in judging the quality of a functional.— 
According to Perdew and Wang,- the exchange potential of a GGA can be calculated as 
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Solid 


u 


J^LDA 








j^RPBE 








j^PBEa=0.52 


Li 


13.0 


15.06 




14.01 




13.51 




13.39 




14.40 


Na 


7.5 


9.19 




7.71 




6.97 




7.48 




8.15 


K 


3.7 


4.49 




3.62 




3.35 




3.42 




3.82 


Al 


79.4 


82.68 




76.49 




74.43 




79.08 




77.90 


C 


443 


467.35 




432.44 




417.49 




449.20 




441.77 


Si 


99.2 


95.13 




87.78 




83.95 




92.74 




90.07 


SiC 


225 


228.70 




211.76 




204.53 




220.72 




216.42 


Ge 


75.8 


72.40 




59.42 




54.07 




67.62 




62.89 


GaAs 


75.6 


74.50 




60.93 




55.44 




69.42 




64.71 


NaCl 


26.6 


32.76 




24.32 




20.56 




24.88 




26.77 


NaF 


51.4 


63.65 




46.73 




40.14 




47.08 




51.12 


LiCl 


35.4 


42.04 




31.81 




26.22 




35.46 




35.20 


LiF 


69.8 


90.66 




71.73 




65.07 




74.12 




76.50 


MgO 


165 


173.41 




147.76 




135.70 




156.54 




154.80 


Cu 


142 


193.22 




143.93 




125.90 




170.78 




156.13 


Rh 


269 


320.43 




263.74 




245.63 




295.16 




276.47 


Pd 


195 


241.40 




180.36 




159.76 




214.94 




194.81 


Ag 


109 


146.51 




96.64 




76.98 




126.06 




109.79 


MRE 




15.27 




-6.17 




-14.84 




0.61 




-0.43 


STD 




13.04 




7.61 




9.52 




8.82 




7.81 


MARE 




16.39 




7.81 




15.27 




6.80 




5.57 




TABLE IV: Bulk moduli Bq (in GPa) 


. MRE, STD and MARE 


as in 


Table |TTI1 




Vx 


dn 


3^'^ 




dF, 
ds 




3 ' 


ds 


dF,' 
ds J 


) (10) 



where u = {2kF)^^n^'^'Vn ■ V| Vn|. The potential thus depends on both the first and second 
derivative of the enhancement form. One advantage of the PBEa form compared to earlier 
attempts to vary the exchange enhancement is that its simpler form leads a smoother ex- 
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FIG. 2: (color online) Derivatives of the exchange enhancement that enter the exchange potential 
(Eq. [10|). Lines for PBEa = 1, WC and mPBE are shown. The thick lines are for the 
term and the thin lines for the 4^ ( | term. 



change potential. This is shown in Fig. [2] where the terms and ^ |^"5~^^J from the 

PBEa = 1, the WC and the mPBE are shown. It can be seen that the WC terms, due 
to the complicated form of X2, have an oscillating behavior while the mPBE converges very 
slowly to zero. The slow convergence of the mPBE is particularly unpleasant as this term 
is multiplied by the 4/3s^, which diverges for large s. Though GGAs in general do not have 
the correct 1/r asymptotic decay, the mPBE does not correct this problem, and the effects 
is an artifact of a poorly constructed functional form. 



A new functional form for GGAs has been proposed, which can be systematically varied 
by the parameter a. One could thereby obtain a functional that combines the virtues of 
semi-empirical GGAs, high precision for a specific set of compounds, with the virtues of 
PBE, obeying the same exact constraints as the PBE. Compared to other semi-empirical 
functionals it is an advantage that only one parameter is used and that the a has clear 
interpretation, which can be used to check whether results are biased towards densely packed 
or open systems. One further advantage of the PBEa is that the simple construction of the 
functional form leads to a smooth exchange potential. 

It is clear, that the limited amount of information in the gradient of the density makes it 





IV. CONCLUSION 
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impossible to make a GGA functional that works equally well for all cases. On higher rungs 
of the DFT ladder,- such as the meta-GGAs, this could be possible. One way to include the 
extra information is make the x dependent on more variables, as in the recently proposed 
TPSS meta-GGA.^^ A different route could be to make a be determined by the density 
and not be the same every where in space. One could e.g. use that for an exponentially 
decreasing density q/p = 1 to identify regions where a should be large. An other idea which 
does not involve higher derivatives was put forward by Wu and Cohen,— who suggested 
making the exchange enhancement dependent directly on the density {Fx{n,s)). 

A subroutine implementing the PBEa can be found at 
http://www.chem.au.dk/~webuorg/georg 



* Electronic address: | geor g@ chem . au . dk 
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